Origins of parton correlations in nucleon and multi-parton collisions 
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We demonstrate that perturbative QCD leads to positive 3D parton-parton correlations inside nu- 
cleon explaining a factor two enhancement of the cross section of multi-parton interactions observed 
at Tevatron at x% > 0.01 as compared to the predictions of the independent parton approximation. 
We also find that though perturbative correlations decrease with x decreasing, the nonperturba- 
tive mechanism kicks in and should generate correlation which, at x below 10 -3 , is comparable in 
magnitude with the perturbative one for x ~ 0.01. 
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Multiple hard parton interactions (MPI) is an impor- 
tant element of the picture of strong interactions at high 
energies. At the LHC energies MPI with p± ~ few GeV 
occur in inelastic collisions with probability of the order 
of one. The issue of MPI is attracting a lot of attention. 

Building up on pioneering works of the early eighties [TJ 
[2], theoretical studies were carried out in the last decade; 
for summary see proceedings of the MPI workshops |3j[4]. 

Multi-parton interactions can serve as a probe for cor- 
relations between partons in the nucleon wave function 
and are crucial for determining the structure of the un- 
derlying event at LHC energies. MPI are currently mod- 
eled by Monte Carlo (MC) generators assuming the pic- 
ture of independent partons. 

A number of experimental studies were performed at 
the Tevatron [5H?]. New measurements are underway 
at the LHC, as MPI may constitute an important back- 
ground for new physics searches. The longstanding puz- 
zling feature of the Tevatron data is that the indepen- 
dent parton approximation which is strongly constrained 
by the HERA data on the hard diffractive vector meson 
production leads to at least a factor of two smaller cross 
section than reported experimentally [8] . To describe the 
data, the MC models employ for the area of transverse 
localization of partons the value which is twice smaller 
than the one indicated by the HERA data. 

Double hard parton scattering in hadron-hadron colli- 
sions contributes to production of four hadron jets with 
large transverse momenta P^^>Aq CD , of two electroweak 
bosons, or "mixed" ensembles comprising three jets and 
7, two jets and W, etc. For the sake of definiteness, we 
will refer to production of four final state jets. 

In [9j [10] we have developed a new formalism to ad- 
dress the problem of multi-parton interactions. QFT de- 
scription of double hard parton collisions calls for intro- 
duction of a new object — the two-particle generalized 
parton distribution, 2GPD. Defined in the momentum 
space, it characterizes non-perturbative two-parton cor- 
relations inside hadron [9]. 

The double hard interaction cross section (and, in par- 



ticular, that of production of two dijets) can be expressed 
through the generalized two-parton distributions 2GPD 

D b a c (x u x 2 ,qlqlA 2 ). 

Here the index a refers to the hadron, indices 6, c to the 
partons, x\ and x 2 are the light-cone fractions of the par- 
ton momenta, and q 2 ,q 2 the corresponding hard scales. 
The two-dimensional vector A is Fourier conjugate to the 
relative transverse distance between the partons b and c 
in the impact parameter plane. 

The two partons can originate from the non- 
perturbative (NP) hadron wave function or, alternatively, 
emerge from perturbative (PT) splitting of a single par- 
ton taken from the hadron. In the first scenario one ex- 
pects that the typical distance between partons is large, 
of the order of the hadron size R, so that the correspond- 
ing correlator in the momentum space is concentrated at 
small NP scale A 2 ~ R~ 2 and falls fast at large momenta 
(exponentially or as a high power of A 2 ). At the same 
time, PT production of the parton pair is concentrated 
at relatively small distances, so that the corresponding 
contribution to 2GPD is practically independent on A 2 
in a broad range up to the hard scale (s) characterizing 
the hard process under consideration. 

Separation of PT and NP contributions is a delicate 
issue. We suppose the existence of a separation scale 
Qo (of order of 1 GeV) such that the A-dependence of 
the correlation function differs substantially for the two 
mechanisms. Once this is done, it is reasonable to repre- 
sent the 2GPD as a sum of two terms: 

D*(x u x 2 ,q*,q*,A 2 ) = 

[2] D b a c (x u x 2 , ql qlA 2 ) + {1] D b a c (x u x 2 , ql ql A 2 ), 

where subscripts and mark the first and the 
second mechanisms, correspondingly: two partons from 
the wave function versus one parton that perturbatively 
splits into two. 

As it has been explained in [10], a double hard inter- 
action of two pairs of partons that both originate from 



perturbative splitting of a single parton from each of the 
colliding hadrons, does not produce back-to-back dijets. 
In fact, such an eventuality corresponds to a one-loop 
correction to the usual 2 — >> 4 jet production process and 
should not be looked upon as multi-parton interaction. 

There are two sources of genuine multi-parton interac- 
tions: four-parton collisions described by the product of 
(PT-evolved) 2 GPDs of NP origin, 

[ 2 ]D a (x 1 ,x 2 ; A) [ 2 pb(x 3 ,X4] -A), (2a) 
and three-parton collisions described by the combination 



[2pa [Ipb + [Ipa [2pb- 



(2b) 



The latter corresponds to an interplay between the NP 
two-parton correlation in one hadron, and the two par- 
tons emerging from a PT parton splitting in another 
hadron — 3 — ?► 4 processes. 

Hard scattering of two pairs of partons in one hadron 
collision event is formally a rare process. Compared with 
production of the same final state in usual two-parton 
interactions (two-to-four), contribution of double hard 
processes (four-to-four) is small. It is suppressed as a 
power of the overall hardness of the process, Aq CD /Q 2 , 
with Q 2 ^pl,M^ ± . 

However, in specific kinematical region the contribu- 
tion of 4 — >• 4 processes may become comparable with 
that of 2 — » 4. This happens in the "back-to-back kine- 
matics" when four final state jets group into two pairs 
each of which has relatively small transverse momentum 
imbalance, £ 2 2± , 5g 4J _ <C Q 2 , with Sij± = Pi± + Pj±- 

In [10] QCD evolution equations for 2GPDS were de- 
rived in the leading collinear approximation, and differ- 
ential distributions in £13, S24 were presented in the form 
resembling the DDT formula for the ^-distributions of 
massive lepton pairs (Drell-Yan process). 

Differential distributions due to 4 — >> 4 and 3—^4 pro- 
cesses exhibit double collinear enhancement: they peak 
at small jet imbalances, S 2 k <C p 2 ± 
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Structure of singularities displayed in Eq. 3a 



(3a) 



(3b) 



m- 

is 



dependent enhancements in two pair imbalances 
typical for 4 — >> 4 processes. The 3^4 processes also 



contribute to Eq. 3a from the region of the transverse mo- 
mentum scales of the splitting n 2 <C max{£ 2 3 , 5 2 ^} ("in- 
ternal splits"). The situation is different when there is no 
QCD emissions between the parton splitting 0^1 + 2 
and the two hard vertices as shown in Fig. [T] This "end- 
point" contribution is enhanced as Eq. |3bJ Singularities 
in Eq.|3]get smeared by double logarithmic Sudakov form 
factors of the partons involved, depending on the ratios 
of proper scales, see [ID] . 





FIG. 1: Origin of 5' 2 singularity in Eq. 3b 



The four-jet production cross section due to double 
parton-parton scattering is conveniently represented as 
a product of cross sections of two independent hard col- 
lisions normalized by the effective correlation area S 
(known in the literature under a misleading name "ef- 
fective cross section"): 



dcr^ da(xi,X2) da(xs,X4) 



1 



(4) 



dt\dt2 dti dt2 S 

It is given by A-integral of the product of 2GPDS Eq. [2] 

I - f d * A 



{[2pa[2pb + [lpa[2pb + [ 2 p a [lpb} • (5) 



The PT 3—^4 end-point contribution to the total back-to- 
back cross section should be added to Eq.|5j see Eqs. (32) 
of [ID]. It does not factorize into the product of two- 
parton distribution functions of colliding hadrons. In 
spite of being numerically small, it is worth trying to ex- 
tract experimentally, as it manifests specific correlation 
between pair jet imbalances, Eq. [3b] 

The NP two-parton distribution [ 2 jD falls fast with A 2 
above lGeV 2 , while depends on A only logarithmi- 
cally. This enhances the contribution to the inverse effec- 
tive correlation area of 3 — » 4 processes (the sum of the 
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FIG. 2: The ratio of 3^4 to 4^4 PT-evolved contributions 
to 2 GPD(xi=x 2 , A = 0) for Q§ = lGeV 2 . 
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second and third term in Eq. [5| by about factor of five 
as compared with the genuine 4 — 4 (the first term) [9]. 

In Fig. [2] the ratio ^D p /^ 2 ]Dp at A = is displayed for 
two gluons with xi=x 2 . It quantifies the strength of lon- 
gitudinal correlation, since the point A = corresponds 
to integral over the transverse distance between partons. 

The numerator was calculated by Eq. 18 of [10] using 
the GRV parametrization of proton pdfs [11]. For the 
denominator we employed the evolution equation Eq. 16 
of [10] by taking for input the product of generalized one- 
parton distributions D(x,q 2 ;t — —A 2 ) according to the 
model of independent partons: 

[2 ]D(a?i, a? 2 ; <Zi > tfi; A ) = D{x u <ft\ A 2 )D(x 2 , g£; A 2 );(6a) 
D(x,Q 2 ,A 2 ) = G(x,Q 2 )-F g (A 2 ). (6b) 

Here G is the standard gluon pdf and F is the two-gluon 
form factor of the nucleon, see [9 j. Fig. [2] indicates that 
the effective correlation area ("effective cross section") 
does depend on the transverse momentum scale, an im- 
portant feature which is not built in in the existing MC 
event generators. 

The following simplified model allows one to get a qual- 
itative estimate of the relative importance of the 3^4 
contribution, as well as to understand its dependence on 
x and the ratio of scales, Q 2 vs. Qq. Imagine that at 
a low resolution scale, Qo? the nucleon consisted of n q 
quarks and n g gluons ( "valence partons" ) with relatively 
large longitudinal momenta, so that triggered partons 
with xi,X2 <C 1 resulted necessarily from PT evolution. 
In the lowest order, a s log(Q 2 /Qo) = the inclusive 
spectrum can be represented as 

G oc (n q C F + n g N c )£, 

where we suppressed ^-dependence as irrelevant. If both 
gluons originate from the same "valence" parton, then 



m D cc ^N c £ ■ G + (n q C 2 F + n g C F N c )e, 
while independent sources give 

[ 2 }D oc n g (n q -l)C F + 2n q n g CFN c + n g (n g -l)N 2 
= G 2 -(n q C% + n g C F N c )Z 2 . 



(7a) 



(7b) 



Recall that the A-dependence is different in Eq. 7a and 



Eq. 7b However, at A = the second terms cancel in the 
sum and we get for the correlator 



D bc (x 1 ,x 2 ;0) 
G b { Xl )G c {x 2 ) 



1 



Nr. 



2{n q C F + n g N c ) ' 



(8) 



The correlation is driven by the gluon cascade — - the 



first term in Eq. 7a It gets diluted when the number of 
independent "valence sources" at the scale Qo increases. 
This happens, obviously, when Xi are taken smaller. On 
the other hand, for large Xi ~ 0.1 and increasing, the 
effective number of more energetic partons in the nucleon 



is about 2 and decreasing, so that the relative importance 
of the 3^4 processes grows. 

It should be remembered that 3^4 mechanism con- 
tributes to the cross section ~ 5 times more than to the 
2 GPD in Eq. ^ Indeed, as Fig. |] shows, at x { ~ 10~ 2 
(Tevatron) the 3—^4 contribution enhances by about fac- 
tor of 2 the four-jet production. 

Thus, an account of 3 — >• 4 processes in combination 
with realistic one-parton GPDs explain the absolute mag- 
nitude of the cross section observed at the Tevatron [SHZ] • 
Fig. [3] is also consistent with the trend observed by DO [7] 
of the increase of 1 / S with p± . A more informative con- 
frontation of predictions with the data would require ad- 
ditional effort on both experimental and theoretical side. 
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FIG. 3: The ratio of 3—^4 to 4—^4 contributions to 4g — >4 jets 
cross section (pa — x) for Tevatron and LHC energies for two 
choices of the starting evolution scale Ql; Q 2 — x 2 s / 4. 



A smooth matching between hard and soft QCD phe- 
nomena at Q 2 ^ 1 GeV 2 allows one to assume that at 
such scale the single parton distributions at small x be- 
low 10 ~ 3 are given by the soft Pomeron exchange. In 
this picture the two soft partons originate from two in- 
dependent "multiperipheral ladders" represented by cut 
Pomerons, see Fig. [4] 
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FIG. 4: 2GPD as two-Pomeron exchange 

The soft Pomeron amplitude is practically pure imag- 
inary. As a result, this amplitude equals that for the 
diffractive cut of the two-Pomeron diagram of Fig. [5] The 
two contributions to the cut are the elastic and diffrac- 
tive intermediate states. The elastic intermediate state 
obviously gives the uncorrelated contribution to 2GPD, 
while the inelastic diffractive cut encodes correlations. 
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FIG. 5: 2iP contribution to 2GPD and Reggeon diagrams 



It is convenient to introduce the ratio 

m A 2x _ [2]D(xi,x 2 ;QlQ 2 ;A) 

[XUX2 ' j " G{ Xl ,Ql)G{x 2 ,Ql) ' 

Then, Fig. [5] allows one to write 

R(xi,x 2 ,t) = R el (x 1 ,x 2 ,t) + R diS (x 1 ,x 2 ,t). 



(9) 



(10) 



The first term (uncorrected partons) is given by the 
product of two single-parton GPDs, 

R e i(x 1 ,X2,t) = gi{xi,t)g 2 (x 2 ,t), (11) 

with gi the ratio of a GPD to the corresponding pdf. The 



second term in Eq. 10 induces parton correlations, both 
transverse and longitudinal. 

By virtue of the QCD factorization theorem for hard 
diffraction [12] , the ratio i?diff / ^ei at t = that quanti- 
fies the strength of the longitudinal correlation is given 
directly by the ratio of the corresponding cross sections. 
The value 

lj = Rdis(xi,x 2 ,0)/ Re\(xi,x 2 , 0) = 0.25 ± 0.05 (12) 

is extracted from the ratio of inelastic and elastic proton 
diffraction cross sections for electro-production of vector 
mesons studied at HERA [15] . For x < 10 -3 this param- 
eter was found to depend only weakly on the incident 
energy (x). Moreover, it stays roughly the same for light 
mesons and J ftp in a wide range of Q 2 , thus confirming 
the hypothesis of a smooth transition between soft and 
hard regimes. 



Determination of the t-dependence of the ratio from 
the data is much more uncertain. Studies of various 
diffractive processes, both "soft" (pp — » p + Mx) and 
"hard" (7 + p -» J/ijj + p, 7* + p -» V + p with 
V = p, lj, 0, J/ip) indicate that the t-dependence of the 
differential cross section is dominated by the elastic ver- 
tex ppJP oc exp(B eX t) with B e} = 5-^6GeV 2 for x < 10" 3 . 

Using an exponential parameterization exp(£?i ne it) for 
the square of the inelastic vertex pMxP, the experimen- 
tally measured ratio of the slopes B- ine \/B e \ ~ 0.28 [T5] 
translates into the absolute value B- lrie \ = 1.4-=- 1.7 GeV 2 . 
The fact that ^-dependence of soft and hard inelastic pro- 
cesses is similar, goes in line, once again, with the logic 
of smooth matching of hard and soft regimes. 

The A integral of the first (4 — >• 4) term in Eq. [5] gives 
the enhancement factor 



(l/5)c 



(1/5), 



l + 2w 



2B e 



inel 



i B e 

B, 



inel 



(13) 



The central value uo — 0.25 yields 77 = 2. Hence, 
the MPI enhancement should persist even for x < 10 -3 
where we expect the PT correlations to diminish. 

In conclusion, we have demonstrated that transverse 
and longitudinal correlations between partons in a nu- 
cleon generated by PT and NP mechanisms play a criti- 
cal role in explaining the observed absolute rate of MPI. 
Dedicated studies of the MPI at the LHC will pro- 
vide a deeper understanding of the nucleon structure 
beyond single-parton distributions. It is pressing now 
for Monte Carlo models of pp collisions to incorporate 
parton-parton correlations and employ realistic single- 
parton transverse space distributions as determined from 
hard exclusive processes at HERA. 
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